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originatea  from  detonation  o;'  1.1  charros  ric'^es  i 

the  center  of  tne  tubes/turnel .  uiessuremcnt  s  of 
pressure  time  history,  employing  piezoelectric 
pressure  transducers,  shoe;  that  front  pro  s- sc  mo  ,  i..,- 
pulse,  positive  duration  ana  time  ol  arrive j.  :cale 
according  to  ooramon  one-dl-r.ensionnl  ,scamj.nf_  ] 
for  small  distances,  L/D,  expressed  in  aiameter 
units,  and  for  longer  distances  gave  neviation:.- 
v;ith  L/D  more  or  less  as  rarar.eter.  The  efforts 
of  vnry-i  nr  wall  rourhne^s  is  .soi,  ini'estmrnted 
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NOMENCLATURE 


C  attenuation  constant  defined  in  text 

D  mo;Jor  inside  tube/tunnel  diameter 

e  averaged  absolute  wall  roughness 

e/D'lOO^  rel.  wall  roughness  in  % 

I  impulse  of  positive  pressure  front 


reduced  impulse 
scaling  factor 

attenuation  constant  defined  in  text 
D( 

thickness  of  plane  charge  in  tube 
distance  charge-observation  point 
side-on  front  pressure  (atm) 
unattenuated  scaled  pressure 

pressure  at  L/D  >  20 


p  (t)  averaged  pressure-time  profile 

cL  w 

Q  charge  weight  (TNT) 

t„  time  of  arrival  of  blast  front 

a 


duration  of  overpressure 


1 


The  success  of  any  physical  investi¬ 
gation  depends  on  the  Judicious 
selection  of  what  is  to  be  observed 
as  of  primary  importance)  combined 
with  a  voluntary  abstraction  of  the 
mind  from  those  features  which, 
however,  attractive  they  may  appear, 
we  are  not  yet  sufficiently  advanced 
in  science  to  investigate  with  pro¬ 
fit. 

J.  Clerk  Maxwell 


I.  INTRODUCTION 

The  need  for  blast  load  prediction  methods  in  tunnel 
configurations  has  long  been  recognized.  This  need 
hac  become  acute  because  of  military  underground 
defence  construction  projects,  civil  tunnel  construc¬ 
tions  for  roads  and  water  and  transport-tunnels  at 
the  factories  for  explosives. 

This  report  presents  the  results  of  a  series  of  ex¬ 
periments  on  one-dimensional  blast  wave  propagation 
in  steeltubes,  diameter  5  cm  and  20  cm,  in  a  con¬ 
crete  tube  of  diameter  80  cm  and  finally  in  a  tiuinel 
in  rock  of  effective  diameter  6  m.  The  blast  waves 
originated  from  detonation  of  TNT  charges  placed  in 
the  center  of  the  tubes/tunnel* 

The  experimental  program  was  devoted  to  a  study  of 
the  validity  of  common  scaling  laws  and  the  attenua¬ 
tion  of  plane  blast  waves.  The  discussion  and  re¬ 
sults  are  presented  in  rather  simple  terms  intended 
to  give  the  reader  a  physical  picture  of  the  phenomena 
which  occur.  The  detailed  discussion  of  the  problem 
would  be  exceedingly  complex  and  beyond  the  scope  of 
this  report.  This  report  is  therefore  not  claimed 
to  be  a  strict  scientific  one,  but  rather  a  survey 
with  some  rules  for  thumb  for  people  in  charge  of 
planning  and  design  of  constructions  mentioned  above. 
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SCALING 
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If  all  parts  of  a  model  have  the  same  shapes  as 
corresponding  parts  of  the  prototype,  the  two 
systems  have  geometrical  similarity.  There  is 
then  a  point-to-point  correspondence  between 
model  and  prototype. 

When  numerical  values  from  measurements  in  model 
are  transformed  to  the  prototype,  we  must  consider 
the  independent  variable  in  the  problem.  A  dimen¬ 
sional  analysis  of  the  relationship  invariably 
leads  to  an  ecuation  of  the  form 

(II. 1)  a  =  f(a^,  .  a^) 

where  the  a's  are  a  complete  set  of  dimensionless 
products. /1/.  In  most  cases  is  it  impossible  to 
get  complete  similarity  in  a  model  test.  There¬ 
fore  we  may  let  some  of  the  independent  variables, 
which  are  thought  to  be  of  secondary  influence 
on  the  phenomenons  to  be  investigated,  diverge  from 
their  correct  value. 

In  two  systems  -  prototype  and  model  -  one  may 
choose  homologous  cartesian  systems  of  axis  (x, 
y,  z)  and  (x' ,  y' ,  z').  Suppose  one  couples 
the  two  systems  together  such  that  homologous 
points  and  times  are  defined  by 

x'  =  k  •  X 

Jv 


Z '  =  ^2  *  ^ 
f  =  kj  •  t 


The  constants  k  ,  k  and  k  are  the  scaling  factors 

Jr  ^ 

for  the  lengths  in  x,  y  and  z  direction.  k^  is 
the  time  scaling  factor. 

The  general  expression  for  similarity  can  be  defin¬ 
ed  from  two  abstract  scalar  functions  f(x,  y,  z, 
t )  and  f'(x',  y*,  z',  t'); 

It  is  a  similarity  between  the  functions  f  and  f 
if  the  ratio  f'/f  is  a  constant  and  the  functions 
are  evaluated  at  homologous  points  and  times.  The 
scaling  factor  for  the  fiinction  f  is; 
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(II. 3)  =  f‘/f 

In  the  evaluation  of  scaling-laws  for  a  system  it 
is  often  useful  to  find: 

a)  The  independent  variables  which  are  of  import¬ 
ance  for  the  experiment. 

b)  Perform  a  dimensional  analysis  of  the  variables 
and  find  invariant  products. 

c)  Pind  the  scaling  factor  for  transformation  from 
model  to  prototype. 

II . 2  The_principle_pf  similarity  in  one-dimensional  blast 
y S  pr9pagation 

The  practical  importance  of  the  principle  of  simi¬ 
larity  lies  in  the  economy  of  effort  it  permits  in 
determining  the  properties  of  blast  waves  and  in 
the  predictions  it  makes  possible  of  the  effect  of 
changed  scale. 

Once  the  form  of  a  function  is  found  experimentally, 
pressure  say,  over  a  range  of  either  distance  or 
charge  size  only,  its  value  for  other  weights  or 
distances  is  known. 

In  order  to  check  the  validity  of  this  proposition, 
let  us  change  the  scales  of  measurements  of  both 
length  and  time  by  a  factor  k^,  i.e.  L2  = 

^0^1’  ^2  *  ^0^1 1*  principle  of  similarity 

asserts  that  the  pressure  and  other  properties  of 
the  blast  wave  will  be  unchanged  if  the  dimensions 
of  the  charge  Q  are  charged  by  the  same  factor  k  , 
or  refering  to  fig  II. 1  the  thickness  of  a  plane® 
charge  in  tube  is  changed  from  1^  to  I2  »  k(,»  l^j. 

This  is  justified  by  the  fact  that  the  Euler  equa¬ 
tions  and  the  Rankine-Hugoniot  equations  describing 
the  blast  wave  if  viscosity  and  heat  conduction 
can  be  neglected,  are  satisfied  by  all  sets  of 
values  of  L,  t  and  1  as  the  scaling  factor  kocancels 
out  as  all  the  derivatives  are  of  first  order. 

This  is  discussed  in  detail  in  ref  /2/, 

Prom  this  we  may  conclude  that  if  the  principle 
of  similarity  is  true  at  any  value  (L^ ,  t^),  e.g., 

pressure  P(r^,  t^)  -  PCk^r^ ,  k^t.])  =  P(r2,  t2),  it 

is  true  for  all  values  of  L  and  t. 
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The  validity  of  the  principle  therefore  depends  on 
whether  similarity  3.s  found  to  hold  true  in  the 
initial  stages  of  the  explosion. 


Fig  II. 1  ]Plane_charges_filling_the_whole_tube 
cross-section 


4 

The  fact  that  the  pressure  and  other  properties  are 
unchanred  if  the  linear  dimensions  of  the  source 
and  scales  of  length  and  time  are  all  changed  in 
the  same  ratio,  does  not  of  course  specify  what  the 
values  are  without  other  informations,  hut  experi¬ 
ments  may  give  the  functional  dependence  The 
principle  can  be  satisfied  only  if  the  pressure 
depends  on  distance  and  time  only  as  a  function  of 
the  ratios  1/L,  t/1  and  hence  the  front pres sure  may 
he  expressed  as 

(11. 4)  p  =  f(i) 

Another  important  property  of  such  a  blast  wave  is 
the  impulse,  For  unit  area  of  the  wave  front  the 
impulse  I  is  given  by 

V 

(11. 5)  I(L,  t')  =  )  p(L,  t)dt 

where  the  origin  of  time  is  taken  to  be  arrival  of 
blast  front  at  L.  The  time  t'  to  which  the  inte¬ 
gration  is  carried  out  should  be  taken  proportional 
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to  the  scaZif\gf  factor  and  we  write  t'  =  A  •  1  where 
A  is  a  function  only  of  1/1^,  and  the  pressure  P 
depends  on  h  and  t  only  hy  the  ratios  lA,  t/1. 

hej'ce 

A(^) 

(II.G)  I  (I,  t*''  =  1  •  J  pC^,  j-)  d  (*) 

*0 

and  trie  integral  being  a  function  only  of  lA»  we 
obtain 

(II. 7)  I  (L,  t)  =  1  •  e;[i,  A(i)] 

where  g  is  an  undetermined  function. 

The  impulse  neasc<red  for  proportional  distance  anl 
'■ime  scales  is  therefore  proper-^ i  o  al  io  f  he  linear 
dimension  1. 

If  one  -^etonates  a  spherical  charge,  Q,  in  a  tube, 
the  detonation  generates  a  blast  v/ave  which  more 
and  more  gets  the  character  of  a  ore-dimensional 
wave  when  it  propagates  dov/n  the  tube  (fig  II.  5)* 


Fig  II. P  Spherical_charges_in_tubes 


Since  the  linear  dimensions  of  a  charge  are  pro¬ 
portional  to  the  cube  root  of  the  volume  and  hence 
weight  Q  the  thieVness  1  in  ^ig  TI.'I  Lr  given  as 
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(11,8)  1®^-^ 


From  the  previous  v/e  are  now  prepared  to  write 
the  scaling  laws  for  the  most  important  blast 
wave  parameters. 


2 

(11,9)  Pressure:  P  *  f 

2 

(II,10)  Impulse:  '"'o  =  K  (^^  ) 


Positive  duration 


(or  duration  of  overpressure): 


(11.11) 


Arrival  time  of  pressure  front  (or  time  from  de¬ 
tonation  to  blast  v/ave  reaches  distance  L) : 

(11.12)  -  1 

P 

Here  f,  g,  h  and  i  are  unknown  function  of  Q/LP  * 
(shortened  equivalent  charge  weight  Qg)  and  may 

hopefully  be  determined  by  eiiperimcnts, 

.i^xperinents  on  this  has  been  done  in  3v/cden 
ref.  /3/. 


II ,  5  i  22_i2_i22S_tH22SiS  • 

A  blast  wave  is  a  physical  discontinuity  which  is 
propagated  through  a  suitable  clastic  fluid  v/ith 
a  characteristic  .Tjpercoriic  velocity  (Mach 
number).  V/hen  the  mcdiuia  exhibits  energy  dissipa¬ 
tive  effects  such  as  viscosity  and  thermal  con¬ 
ductivity,  a  reduction  in  the  magnitude  of  these 
discontinuities  and  an  increase  in  the  length  of 
the  discontinuous  region  occurs. 

An  added  attenuation  is  introduced  when  the  blast 
wave  is  confined  to  the  interior  of  a  bounded 
medium  such  as  a  tube.  The  added  factor  nay  be 
separated  into  tv/o  contributions,  one  being 
characteristic  of  the  fact  that  the  blast  wave 
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Is  generated  froe  detonation  of  more  or  less 
■pherlcal  charges  within  the  tube  Itself,  and 
the  other  being  characteristic  of  the  walls  and 
the  strength  of  the  blast  wave. 

The  first  type  of  attenuation  gives  a  pressure  de 
crement  which  Is  very  steep  and  is  in  effect  the 
first  few  tube  diameters  in  traversed  distance 

(fig  II. 5). 


Fig.  II. 3*  Blast  wave  formation  inside  a  tube 


The  pressure  across  the  shock  front  attenuates  as 
it  travels  down  the  tube  due  also  to  the  frictional 
losses  on  the  walla  of  the  tube.  For  very  short 
tunnel  lengths,  i.e.  for  tunnels  whose  lengths  are 
not  very  many  times  greater  than  the  diameter  '.oi 
its  average  cross-sectional  diameter  for  the  tunnel) 
attenuation  due  to  wall  roughness  is  very  small. 
However,  for  long  tunnels  having  lengths  equal  to 
about  50  or  more  diameters,  significant  attenuation 
of  the  blast  wave  should  be  expected.  One  of  the 
primary  aims  of  this  program  is  to  experimentally 
determine  the  magnitude  of  this  attenuation. 


Fig,  II. 4.  • 
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Without  wall  attenuation  and  different  other  dissipf.tive 
effects  for  different  tubes/tunnels  and  charges,  the  ex¬ 
perimental  points  for  front  pressure  versus  equivalent 
charge  after  eq  (II. 9)  should  lie  on  one  smooth  curve  as 
shown  in  fig.  II. 4.  Attenuation  due  to  wall  roughness 
would  move  this  curve  with  L/D  more  or  less  as  parameter 
for  one  type  of  tube  and  roughness.  This  comes  about 
from  principle  of  similarity  described  in  part  11,2, 
which  evidently  fails  in  any  circumstances  for  forces 
not  scaling  geometrically.  One  such  force  is  the  effects 
of  viscosity  which  gives  rise  to  terms  of  second  order 
derivatives  in  the  equation  of  motion,  and  with  this  the 
equation  of  motion  is  not  satisfied  on  substituting 
P(koLl,koti)  etc  for  p(Li,t-))  etc. 

At  the  last  front,  processes  similar  in  effect  to  mac¬ 
roscopic  viscosity  and  thermal  conduction  must  be  signi¬ 
ficant,  but  this  turns  out  to  be  of  no  practical  signi¬ 
ficance  to  the  applicability  of  the  principle  of  simi¬ 
larity. 

For  the  tunnel  the  effective  diameter  is  taken  as  D  = 
(4A/lI)i,  where  A  is  the  tunnel  area,  applicable  for 
scaling  of  length  in  frictionless  flow.  Friction  effects, 
however,  depend  on  the  length  of  the  tunnel  periphery  C 
in  contact  with  the  moving  gas,  and  the  effective  dia¬ 
meter  is  more  properly  C/II.  In  our  case,  however,  the 
tunnel  diameter  was  not  defind  better  than  difference 
between  the  tv/o  definitions  because  of  somewhat  varying 
tunnel  diameter. 


Ill 
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TEST  PROGRAM 

The  tubes  used  were  chosen  as  to  make  a  rough 
test  of  the  validity  of  common  scaling  laws.  As 
our  Office  of  Test  and  Development  neither  had 
the  money  nor  the  time  to  make  extensive  tests, 
only  steeltubes  of  diameter  5  and  20  cm,  a  con¬ 
crete  tube  of  diameter  80  cm  and  a  tunnel  in  rock 
of  effective  diameter  ^00  cm  with  constant  rough- 
i^ess  for  each  tube/tunnel  were  used. 

General  test  descriptions  for  these  series  of 
experiments  are  tabulated  in  table  III.1. 


Tube/Tunnel 

Steel 

Steel 

Concrete  Rock 

DiSideter 

D(cm) 

5 

20 

80  600 

Wall  rough¬ 
ness  e/D  100% 

9,5» 

(e<'0,26  hd) 

0,13%  0,03%  8% 

(e'*0,26  inm)(e«0,2  inm)(e*0,5  m) 

Charges  TNT  \ 

(g)  between*' 

1,5-17,5 

1,5-200 

9,5-5000  21,5-92300 

Distance  charge 
-  seas.  L/D 
between 

6,4-705,2 

4,5-150,7 

4,85-100,5  2,92-38,33 

Side-on  over¬ 
pressure  p(ato) 
between 

0,2-77 

0,2-32 

0,07-23  0,005-1,5 

x)  Electrical  detonation  cap  no  6,  equivalent  TNT  charge 
1,5  g  included. 


Table  III.1 


Test  description 
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III. 2 


III.  5 


Pencil  shaped  Lead  Zirconate  transducers  were  used 
in  conjuction  with  an  oscilloscope  and  polaroid 
camera  to  measure  pressure-time  history.  These 
transducers  were  centermounted  in  the  tubes/tunnel. 

Pig  II. 5  shows  the  80  cm  diameter  concrete  tube 
andll.6  a  portion  of  the  3  cm  diameter  steel  tube. 


Rpugtoegs  of  _tubes_and_  ti^el  in  rock 

Tunnel  wall  roughness  was  estimated  by  taking  the 
percentage  of  the  ratio  e/D,  where  e  is  the  roughness 
particle  diameter  and  D  is  the  tube/tunnel  dia¬ 
meter* 


e  (mean  value) 


Pig. III. 1.  Definition  of  roughness  ^  •  100% 


As  for  the  steeltubes  the  roughness  originated 
from  corrosion.  Por  the  concrete  tube  the  rough¬ 
ness  was  the  standard  production  roughness  for 
PREMO-tubes,  manufactured  at  A/S  Dalen  Portland- 
Cementfabrik,  Brevlk,  Norway. 

The  tunnel  used  was  an  abandoned  railway  tunnel 
with  raw-detonated  rock.  The  tiinnel  had  a  curva¬ 
ture  of  radius  approximately  100  m.  It  is  thought 
that  the  effect  of  this  bend  can  be  neglected  and 
the  tunnel  may  be  considered  as  being  straight. 


55El25ll2_S^2rB22 

The  explosive  used  were  cylindrical  pressed  TNT 
charges  with  a  hole  drilled  along  the  cylinder- 
axis  to  fit  the  electrical  detonation  cap  no  8. 
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Several  charges 


taped  together 


Fig  III. 2. 


In  a  aeries  of  calibration  testa  the  electrical 
detonation  cap  was  found  to  he  1,5  g  equivalent 
TNT  charge  as  regarde  front  pressure. 

The  charges  were  hanging  in  center  of  tube  in  a 
string  going  through  a  small  hole  on  the  upper 
side  of  the  tube.  In  the  tunnel  teats  the  heavier 
charges  were  placed  on  a  wooden  platform  as  to  be 
placed  roughly  in  the  center  of  the  tunnel.  The 
different  shaped  charges  is  not  believed  to  have 
had  any  significant  influence  on  blast  waves  at 
longer  distances  5). 


Ill, 4  Transducers 

The  gages  used  were  IC-35,  from  Atlantic  Re¬ 

search  Corporation.  This  is  a  pencil  shaped  Lead 
Zirconate  and  side-on  pressure  measuring  type 
transducer.  The  nominal  characteristics  in  table 
III, 2,  Typical  pressure  time  profile  from  test 
in  shock  tube  is  shown  in  fig. III. 4. 

The  transducers  were  calibrated  at  the  factory  to 
within  *  1%  reproducibility  in  the  pressure  range 
from  1,7  ato  to  15.6  ato. 

A  best  fit  calibration  curve  was  used  by  the  com- 

Suter  in  evaluation  of  frontpressure  and  impulse, 
p  to  5%  uncertainty  in  calibration  was  deemed 
adequate  in  view  of  the  larger  errors  introduced 
in  interpreting  the  oscillographs. 


Model  LC-33 


Fig  III. 3.  ^i262sions_of_LC-33_tr^sducer 


jable  III. 2 

fv'OMirJAl  CHAKACTFRISTICS 

Open  Circuit  Receiving  Sensitivity  (1  Keps) 

db  reference  1  volt  microbar .  -101 

,  Volts;  psi .  0.6 

Charge  Sensitivity  (micromicrocoulombs  psi) .  2500 

Frequency. Response,  flat  i-2  db . 1  cps  to  80  Keps 

Capacitance  (microfarads)  0.004 

DC  Resistance  (megohms) .  >500 

Maximum  Static  Pressure  (psi) .  >500 

Maximum  Transient  Pressure  (psi) .  ~1000 

Temperature  Range  (“C) .  -40  to  +100 

I  meant/ 

Pressure  *2%  (,  s;1 .  100 

Temperature  +3%  {°C) . .  -60  to  +70 


ORIENTATION . FACE  ON 

MEDIUM . NITROGEN  GAS 

PRESSURE . 38  PSI 

VERTICAL . 5.0  V/div 

HORIZONTAI _ 0.02  ms/div 


Fig  III. 4  TYPICAL  PRESSURE-TIME  PROFILE 
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III.  5  ?£§5§^licer_mo\mting 

The  way  in  which  the  gage  was  mounted  in  the  5 
and  20  cm  tube  is  shown  in  fig.  III. 5 


Fig .  Ill .  5 .  Tr^sducer_moi^t  ed_in_^_and_20_cm_tube 


For  the  5  cm  and  20  cm  tube  only  one  transducer 
v;as  mounted  in  the  same  tube,  but  as  for  the  80  cm 
tube  four  were  mounted  at  the  same  time,  fig. III. 6. 

The  transducer  body  and  mounting  tube  are  of 
appreciable  diameter  compared  with  that  of  the 
smallest  tube  tested.  The  area  constriction  due 
to  the  3/4-in.  mounting  tube  in  this  case  is  about 
14  percent.  The  presence  of  a  gage  can  locally 
affect  the  peak  pressure  and  impulse  under  these 
conditions. 


3/4"  tub* 

Fig  III, 6.  Tran8ducer_mo^ted_in_80_C5_tube 

In  the  ttinnel  5  transducers  were  mounted  at  the 
same  time  as  shown  in  fig.III.7« 


Fig  III, 7,  Transducer  50^ted_in_6_m_t^nel_in_rock 


III ,  6  rument  at  ion 

The  electrical  arrangement  shown  in  fig  HI. 8  was 
followed  for  the  5  and  20  cm  tubes.  The  Kistler 
Charge  Amplifier  v Model  566)  served  in  most  cases 
to  match  the  high  impedance  of  the  transducer 
(  >  10®  ohm)  to  the  oscilloscope,  a  Tektronix  Model 
502  A,  A  Polaroid  type  oscilloscope  camera  recorded 
the  trace.  In  many  cases  a  large  capacity  served 
as  a  timeconstant  increasing  device. 

Triggerli*  was  done  by  taping  a  loop  of  insulated 
wire  around  the  charge.  In  detonating  the  charge, 
a  well-defined  break  of  wire  is  obtained.  The 
uncertainty  intriggeringis  estimated  to  less  uhen 
500  /A5  . 


Fig  III. 8.  §chematic_electrical_flow_chart^for_5 
and_^20^cm_diameter_tube_test8 

For  the  less  transient  blast  waves  in  the  80  cm  tube 
and  6  m  tunnel',  a  CEO  recording  oscillograph  (Con¬ 
solidated  Electrodynamics)  was  used  with  preferably 
fluid  damped  galvanometers  type  7-565,  fig  III. 9. 

The  galvanometer  has  the  following  specifications: 

Undamped  natural  frequency  1670  Hz 

Flat  (5^)  frequency  range  0-1000  Hz 

Undamped  DC  sensitivity  2  mA/inch 

Terminal  resistance  69  ohm 
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Transducers 


f  ? 

Kistler 

CEC  -  5  -  "12^ 

oscillograph 

7  -  565  type 
galvanometers 

5  channels 

80  cm  tube 

7  channels 

6  m  tunnel 

Triggering 

CUclPgU  ■'  ' 

amp  1. 566 

Firing 


A Start /St op 


Control  unit 
syncronmotor 
with  micro- 
switches 


Pig  III, 9.  Schematic_electrical_flow_chart_for 


in  rock  tests 


Typical  records  from  oscilloscope  and  oscillograph 
is  shown  in  the  reproductions  in  fig  III. 10. 


/ 

♦ 

T 

Q  =  9. 


c’",  L  -  ''^2^2  cr 


3/div,  Y  =  0,3^.  ato/div 


1  / 

+  '  .  1 

I 

’1  i  /  ■ 

1  S 

- 

. f. . 

.  .  'M  , 

:  ''1  . 

. . 

•  s 

\A.Jk.A.AA-^ 

1 

20  cm,  L  =  5C7  cm 

25,5  G  TiYT 

1  iTiG/div,  Y  -  0,67  ato/di 


•I  I  I  I  I  I  ^Tvw^ 

D  -  60  cm,  L  =  44,85  -m,  Q  =  500  g  TNT 
X  =  100  m.G/div,  Y  =  0,068  ato/div 


III  I  I  I  .  r 

D  =  o  m,  L  =  27  m,  Q  =  7  Ng  TNT 
X  =  100  ms/div,  Y  =  0,23  ato/div 


Fig. Ill . 10  T22ical_records_fro2_oscillosco2e_^d_oscillogra2h 
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aurements 

As  frontpressure  is  uniquely  determined  from  tabu 
lated  blast  velocity  data,  a  simple  check  of 
frontpressure  measurements  from  pressure-time 
records  is  possible  by  simple  velocity  measure¬ 
ments.  This  was  done  with  the  5  cm  diameter  tube 
5  mm  diameter  holes  were  drilled  through  the  wall 
with  approximately  20  cm  interval  and  selfmade 
folioswitches  were  taped  over  the  holes  and  con¬ 
tact-time  interval  was  recorded  by  a  Hewlet- 
Packard  coxinter.  Comparison  of  the  two  methods 
is  done  in  fig.  III. 11,  and  the  difference  is 
not  more  than  expected  from  the  time  of  arrival 
measurements. 
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Fig.  Ill  .11 
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IV.  PREPARATION  OP  DATA 

IV, 1  inter2retation_of_trace8 

As  can  be  seen  in  fig  III. 10  the  records  posessed 
an  irregular  pressure-time  trace,  because  of  noise, 
irregularities  in  the  blast  wave  because  of  non¬ 
ideal  one-dimensional  flow, overshoot  in  transducer, 
oscillograph  etc.  The  records  were  then  "averaged" 
by  drawing  a  best  fit  smooth  curve.  However,  the 
determination  of  pulse  height  by  this  method  is 
estimated  to  give  at  least  uncertainty.  For  the 
determination  of  the  impulse  the  uncertainty  is 
even  more,  approx.  20%  partly  because  the  positive 
duration,  t^,  is  very  \incertain.  Definitions  are 

given  in  fig  IV. 1, 


P 


dt 


t 


Pig  IV. 1.  5eTinition_of_garameter8 


The  method  for  conversion  of  the  pressure-time 
traces  to  paper-tape  is  described  in  detail  in 
ref  /6/.  The  technique  is  quite  simple  and  is 
shown  schematically  in  fig.  IV. 2. 
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Picture 


Pig  17.2.  §2bcmatic_conversion_of_2i2tures_to 
gager-tape 

Two  potentiometers  give  voltages  which  is  a  measure 
of  the  X-  and  y-coordinates  of  a  sufficient  chosen 
number  of  discrete  points  on  the  averaged  trace 
(typical  20  points  were  chosen). 


The  further  modification  of  converted  data,  main  pro¬ 
grams  etc.  are  too  lengthy  to  be  described  in  detail 
here,  but  some  features  will  be  mentioned.  The  com¬ 
puting  was  done  with  an  UNIVAC  1107. 

The  frontpressure  P  was  picked  out  as  the  largest 

ordinate. 

The  impulse  of  positive  pressure  front: 
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t^(p=0) 


I 


Pav  ‘ 


was  computed  with  common  trapezoid-integration. 


t_^  was  found  as  time  from  pressurefront  to  averaged 
pressure  curve  again  reaches  zero. 


V. 

V.1 
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RESULTS  AND  DISCUSSION 

Fig  V.1  to  V,5  presents  the  data  relating  the 
observed  front  overpressure  p  to  tube/tunnel 
position  in  L/D  units  with  TNT  charge  weight  as 
parameter. 

From  this  is  it  clear  that  for  the  same  tube,  the 
higher  the  peak  pressure,  the  greater  the  blast- 
wave  attenuation. 

At  longer  b/D  distances  the  results  demonstrate 
that  the  attenuation  is  described  roughly  by 

k^(p„) 

(v.1)  p  .  C(pg)(L/D) 


Taking  the  initial  pressure  p^  at  L/D=20,  we  may 

find  the  constants  k(p2Q)  from  the  averaged  curves 

in  fig  V,5,  (by  extrapolating  straight  lines  back¬ 
wards  to  L/D«20  for  the  5  cm  tube) . 

This  is  shown  in  fig  V,6.  For  p^q  ^  'I  atm  the 

functional  dependence  of  k^j  is  described  by 


(V.2)  k^  at  "0,455  log  P2q  -  0,655 

and 

(V.5)  C  ^  6,79  . 

hence 

fv  a')  T.  70  ,L  log  P20  -  0,655 

(V.4;  p  =  6,79  P2q  ’ 


This  is  shown  in  fig  V,6  and  V.7  to  be  rough 
approximations  for  the  5»  20  spd  80  cm  diameter 
tubes  for  P2q  in  the  region  1  t  o  200  atm  (extra¬ 
polated  back  from  "long  range  attenuation"  whereas 
for  the  attenuation  tunnel  in  rock  and 

for  the  80  cm  diameter  tube  we  do  not  have  this 
simple  functional  dependence. 

The  exponent  k  seems  to  reach  a  more  or  less 
''tent  value  for  P20  ^  atm. 
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As  stated  before,  it  is  evident  that  attenuation 
is  a  function  of  the  relative  wall  roughness, 

i.e.  dependent  on  the  parameter  ^  .  Within  the 

experimental  uncertainty  one  would  not  expect 
significant  difference  in  attenuation  for  the 

steeltubes  and  the  concrete  tube.  between 

0,03%  and  0,5%)«  This  is  confirmed  in  fig  V.6. 

A 

The  tunnel  in  rock  with  roughness  jj  8%  seem  to 

give  more  attenuation  in  the  pressure  region  of 
overlapping,  i.e.  V20~  ^  atm. 

Additional  experiments  with  systematic  change  of 
rougliness  in  one  tube  and  comparison  of  different 
tubes  with  varying  diameter  are  hoped  to  be  the 
next  step  in  our  series  of  experiments. 
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Fig.V. 1  Frontgressure_versus_tube_£Osition_in_L/D_units 
in  ^ 
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n=20cm 


1  >  408  10*  2  468  10*  2  468  10* 
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Fig .  V.  2  Front2ressure_versus_tube_£osit  ion_in_L/D_^it  s 
in  20  cm  tube 
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Fig .  V .  5  l!r2S^Er®  ® 

in  80  cm  tube 


-  29  - 


mwMi 


ssssx8;;a;i;«as  win  ia:;3giinnttM«iii|in  ins  same  n»8a;imBaatti»ttnin«aeannanau>MnBn 
ssi:xmsai8B::»nns:asiiiiHnBsntBunra!ss:i»nm»Bcniinnanni»in8»ia8»BsaBBBBi 
siwmcaBcaganaaaMKdBititittWttttiaiiiBBagaHMPBWiBBiaattBWBaasBBBBi 
ssr^TaaggBnaasntiBnnuKWBsaBaeangBBBBiiiiiiniiiiiriininaBBBPBBtt 


iBBrnwinawa,.  aiBfcJi^wBfc*iwi».^Ba»a;Bai 
aBnsmBBsnttmisBss  ?inPB?mnHBHinBBaai 
i8ngtta:gB;gim»s«C8gngBBi>3gB«ik*wgBgi 
iiaaa8i8C’CB«r>gBBt;!<HigB»iMBi 


liiiimilllli 


8  10* 

3 

!  4 

1  0  1 

1  10*  2 

Fig.V.4  Frontgressure_versus_tunnel_£Osition_in  L/D  units 
in  6  m  tunnel  in  rock 
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Checkinp;  of  the  validity  of  the  scaling  law  proposed 
in  chapter  II.?  eq  (II. 9)  for  the  front  pressure 

(II. 9)  p  .  r(-^)  .  f(Qv) 

LD^  ^ 

is  shown  in  figures  V.P  to  V.12. 

The  5  cm  tube  resul‘»'.'=  in  fip  V.8,  20  cm  in  fig 
80  cm  in  fie  V.10  and  the  6  m  tunnel  in  rock  in 
fig  V.'''"!.  The  influence  of  the  wal]  attenuation 
clearly  manifests  itself  in  the  lowering  of  the 
pressure  curves  with  increasing  L/D  values.  Direct 
comparison  of  all  data  is  given  in  fig  V.'l2. 
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The  overall  feature  is  that  the  deviation  from  the 
ideally  expected  scaled  sniooth  frontpressure  curve 
is  greater,  the  smaller  the  tube  cross  section. 

Plotting  the  pressure  versus  L/D  v.'ith  as  para¬ 
meter  we  may  extrapolate  to  L/D  — ►  0  and  get  the 
"unattenuated"  scaled  pressure  curve,  PqqCQj)* 

The  average  L/D  influence  on  the  attenuation  is 
evaluated  to  have  exponential  behaviour  as 

(V.6)  p  =  Pg^exp  (-k  •  (L/D)) 

The  damping  constant  k  depends  on  p  and  tube/ 

tunnel  and  is  expressed  as  the  fractional  loss  in 
frontpressure  per  diameter  of  travel,  i.e.,  as 


(V.7)  -D(dp/dx)/p  =  k 


This  constant  is  estimated  in  fig  V.15. 

The  tunnel  shows  considerably  greater  wall  attenua¬ 
tion  than  the  smoother  tubes. 

Enerich  et  al.  in  ref  /4/  gives  a  value 

k  =  (1,4  i  0,8)  •  10“'^  as  a  mean  value,  for  a  5,5 

and  15,5  cm  diameter  tube  with  p„  between  0,1  and 
’  so 

1  atm,  which  roughly  lies  in  the  region  of  our 

estimations. 


40 


rig  V.I5  EE22^i25Sl_fl?:?fi§’®®_i2_£S®ssure_froo  scaled 

Tjmattenuated  2Hrve  per  xmit  distancg  of  travel 
ln_tubf;_diameter£ 

Comparative  tests  done  by  Norwegian  Defence 
Research  Establishment  /5/.in  a  tunnel  in  rock 
of  effective  diameter  4,8  -  0,5  m»  i.e.  rough¬ 
ness  about  15%  with  charges  ranging  from  10  to 
500  kg  TNT  and  L/D  =  1,9  -  12,  shows  good 
agreement  with  our  experiments. 


V . 2  ISEHiS® 

The  scaling  law  proposed  in  chapter  11,2  eq. 
(11.10)  for  the  impulse  1 
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(11.10) 


hereafter  called  reduced  impulse) 


Is  checked  with  experiments  in  fig  V.14  -  V.18. 

The  data  here  are  considerably  fewer  than  for  the 
pressure  which  is  because  rarefaction  waves  from 
the  open  end  of  the  tubes  sometimes  destroyed  the 
correct  pressure  time  history. 

Direct  comparison  at  comparable  distances  in  L/D 
units  in  fig  V.17  and  V.18  shows  a  significant 
lower  reduced  impulse  for  the  tunnel  in  rock  at 
large  L/D  (“  38)  compared  with  the  tubes. 

We  may  use  the  same  procedure  as  described  in  V.1 
for  the  pressure  letting  L/D  — ►o  and  find  the 
"unattenuated ”  scaled  impulse  cui’ve.  For  L/D  ^  7 
the  attenuation  is  not  significantly  different  for 
the  tubes  and  tunnel  in  rock. 
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V.5 

The  scaling  law  proposed  in  chapter  II. 2  eq 
(II. 11)  for  the  positive  duration  t^ 

(II.11)  ■'*'  Q -  «  ^ 

is  checked  with  experiments  in  fig  V.19. 

We  are  seeing  the  same  tendency  here  as  for  the 
impulse  -  lowering  of  the  reduced  positive  dura¬ 
tion  curves  with  increasing  L/D  values.  The 
wall  roughness  seems  not  to  have  a  simple 
scaled  influence. 

For  I/D  ^  7  the  positive  duration  seems  to  scale 
roughly  for  all  tubes  and  tunnel. 

We  may  propose  here  too  that  loading  analysis 
be  made  utilizing  the  I,/D  ^  7  data  only  to  have 
a  safety  factor. 
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V.4  ?iS®_2£-S££iYii 

Time  from  detonation  to  blast  wave  reaches  distance 
L,  was  found  to  scale  as 


(11.12) 


i 


which  is  checked  with  the  experiments  in  fig  V.20 
for  all  tubes  and  tunnel  in  rock.  We  have  no  signi¬ 
ficant  deviation  in  this  case  from  one  smooth  curve 
with  increased  L/D  values  witnin  the  experimental 
uncertainty 


For 


and 


10“^;^  < 
~  1  r,2  ~ 


the  scaled  curve  may  roughly  oe  described  by  the 
straight  line  in  the  log-log  plot  as: 


(V.8) 


-1,05 


or 


jj0,1  .  ,1,05 
(V.9)  t,^(ms)  c-  2,0^^ - - 
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VI  SUMMARY 

Frontpressure  of  an  one-dimensional  blast  wave 
generated  by  a  charge  detonating  in  a  tunnel 
scales  as 

(VI. 1)  p  =  f  (-%) 

with  an  attenuation  dependent  on  wall  roughness 
and  L/D  parameter.  The  average  L/D  influence 
on  the  "attenuation”  from  an  "unattenuated" 
scaled  pressure  curve »  ^so^^E^  (VD — >0),  is 

evaluated  to  be 


(V.2)  P  =  Pg^exp  [-k(L/D)] 

Estimates  of  k  is  given  in  fig  V.15  for  different 
tubes  and  tunnel  in  rock. 

Typical  for  1  atm,  20  and  80  cm  tubes 

S  w 

with  wall  roughness  0,03  -  0,3%  have 
k  a  (4  i  2)  •  10“^. 

For  the  6  m  diameter  tunnel  in  rock,  roughness 
^8%,  k  ^  (3  -  1)  •  10”^. 

The  reduced  impulse  scales  as 

(VI. 5)  =  g  (-%) 

with  an  L/D  influence  dependent  on  roughness 
similar  to  frontpressure. 

The  reduced  positive  duration  scales  as 
(VI. 4)  ^4-  • 

5 


h  (~^) 
LD^ 


v/ith  similar  effects  as  described  above. 
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The  reduced  time  of  arrival  scales  as 

t  .  Q 

(VI. 5)  — '-r -  =  i  ^“2^ 

which  does  not  seem  to  have  a  significant  L/S 
dependence 


For 


tg^(ms)  —  2,05 


j50,1  .  j^1,05 

- ^ - 


with  D(m) 

L(m) 

Q(kg  TNT) 


It  is  recommended  that  loading  analysis  to  be  made 
utilizing  the  “smooth  tube”  or  lim  L/D— >0  data 
only.  It  is  felt  that,  although  roughness  does 
not  seem  to  scale,  any  roughness  in  a  tunnel  would 
tend  to  be  a  safety  factor  If  doors,  equipment  etc 
are  designed  to  withstand  “smooth  wall”  loads. 

The  "smooth  wall”  reduced  parameters  for  the  blast 
wave  are  presented  in  fig  71.1  with  rough  estimates 
of  the  uncertainty. 
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